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Abstract 

Recently observed c-axis optical sum rule violations indicate non-Fermi liq- 
uid in-plane behavior. For coherent c-axis coupling, the observed flat, nearly 
frequency independent c-axis conductivity cri{uj) implies a large in-plane scat- 
tering rate F around (0, vr) and therefore any pseudogap that might form at 
low frequency in the normal state will be smeared. On the other hand inco- 
herent c-axis coupling places no restriction on the value of F and gives a more 
consistent picture of the observed sum rule violation which, we find in some 

cases, can be less than half. 
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C-axis electrodynamics is important in distinguishing high-Tc cuprates from conventional 
superconductors and for understanding mechanism in the cuprates Since the cuprates 
have layered structures, interlayer coupling between two adjacent Cu02 planes plays an 
essential role, and its exact nature impacts on the recent experimental observation of a 
significant violation of the conventional optical sum rule of Ferrel, Glove and Tinkham . 
The normalized missing spectral weight (NMSW) is different from one and of order of a half 
in several high-Tc cuprates. NMSW is the difference between the area under the real part 



of the conductivity cr"'^'^-' (cj) = ai^^\uj) + ia2^^\uj) in the normal minus superconducting 
state, devided by the superfiuid density which is obtained from aKtu). We have pointed 
out in Ref. ^ that in order to explain NMSW of less than one, it is necessary to consider 
non- Fermi liquid models for the in-plane dynamics, regardless of the nature of the interlayer 
coupling. A possible such model introduces a pseudogap in the normal state above Tc as 
is observed 0. Another is the "mode" coupling model determined from consideration of 
ARPES data by Norman a/ 10. This model has been used to describe kinetic, as opposed 
to potential energy driven superconductivity. 

In this paper, we show that coherent c-axis coupling, even with a pseudogap, cannot 
easily explain recent experimental findings, but that incoherent coupling describes them 
including the observed low frequency behavior of the effective NMSW. 
The interlayer coupling is represented by a Hamiltonian 



Hc = Yl [tij4i^Ci2a + H.c] , (1) 

where the hopping matrix tij describes weak interlayer tunneling and c^^ creates an electron 
with spin | at the site i in the plane 1. We classify interlayer couplings in two classes because 
the results are remarkably different depending on their nature: i) One is coherent coupling 
{tij = tj_), which originates from an overlap of electronic wave functions between the two 
planes. For an in-plane Fermi liquid, it was shown in Ref. M that the superfiuid density (p^) 



is equal to the missing spectral weight (A^„(c<Jc) —Ns{uJc) = 8 Jq+ duj' (t"^(co'') —al^{uj' 



where 



ujc is a cutoff frequency of the order of a bandwidth). Thus coherent coupling (unless the 



density of states has a strong variation with energy) can explain the experimental results on 
optimally doped YBCO and over-doped T12201 |]10[. ii) The other is incoherent coupling, 
for which tij = ViSij, where Vi is an impurity potential. In this case, we showed in Ref. 0] 
that NMSW > 1.58. The characteristic difference between coherent and incoherent coupling 
is whether or not electron momentum is conserved in the interlayer transfer. 

In the presence of an external vector potential Az, He is modified to Hc{Az) by the phase 
factor exp{±ieAz). It is sufficient to expand Hc{Az) to second order in Az- The current 
jc = -6Hc{Az)/SAz = jp + jd, where jp = -ied tj- [ctiaCi2a - ct2aCiia] and jd = e^d?HcAz 
with d the interlayer spacing. In linear response theory, (jc) = [— II + e^(i^(ifc)]^z, where 11 
is the current-current correlation function associated with jp and {H^ is the perturbation 
of id due to He- 

From the c-axis conductivity sum rule P,pJ]| the superfiuid density ps can be written as 



p, = 8 rrf^[<H-aJ,Hl -47reV[(ifc)'-TOi, (2) 

where uj^ is the cutoff frequency for the interband transitions that He does not describe, and 
we use units such that h = c = kB = I and set the volume of the system to be unity. 

The penetration depth Ac can be calculated in two ways. Based on the Kramers-Kronig 
relation for the conductivity, we obtain Ac as l/47rA^ = limj^_^o[^lKicrc(0, cu)]. Alternatively, 
using Eq.(0) we can also calculate Ac(= l/i/p7). Combining these two equations for Ac, we 
obtain the formula: 

{Nn-Ns) _l l Eu.J:^.,p\tk-p?[Gik,co)G{p,cu)-Go{k,u)Goip,uJ)] 
Ps 2 2 E.Ek,p|ik-pPF(k,a7)F+(p,^) ' 

where G(k, cu) and F(k, cu) are superconducting Green functions and G'o(k, cu) is in the 
normal state. For coherent c-axis coupling (|tk~pP = ^l^k-p) electron momentum parallel 
to the plane is conserved and may depend on the in-plane momentum. For incoherent 
coupling, an impurity configuration average is imphed over a potential Vi (|tk-pP = iHi-pP), 
and electron momentum is not conserved, and k and p remain unconstraint. 

Coherent coupling. As the simplest case, one can consider a normal state spectral function 
Aq{]s.,uj) = T / [{uj — ,^k)^ + r^] based on Fermi-liquid theory. The in-plane scattering rate 



r is expected to depend strongly on direction of the momentum k in the two dimensional 
Brillouin zone. Cold spot exists along (vr, vr) and hot spot along (0, vr) p. The coherent c-axis 
matrix element t±, however, is itself dependent on direction and t\ varies as cos^(2(/)), where 
(f) is the angle defining the in-plane direction of k. This factor makes the c-axis conductivity 
sensitive mainly to the hot spot, although it remains Drude-like. Here we can ignore both 
dependences and interprete F as the scattering rate coming from (0, vr) region (hot spot). 
However, except for YBCO at optimum doping, ai{uj) vs u is found to be almost flat over 
an energy range of a few 100 meV. This implies that F and F^ need to be several 100 meV. 
Non-Fermi liquid models with large F have appeared in the literature based on ARPES data 
in Bi2212. Norman et al. have derived a simple phenomenological model - called the 
mode model in which the normal state is described by a large frequency independent F and 
the superconducting state by T{uj) = F 9{uj — A(0) — Umode), i-G. a Fermi liquid is recovered 
in this state as T^u) = for < A(0) + ujmode, where A(0) is the gap and Umode is the 
frequency of some collective mode. In addition, a pseudogap A can be introduced in the 
normal state as in the preformed pair model []T2[ with no long range phase coherence above 



Tc. We take the pseudogap to have the same d-wave symmetry as does the superconducting 
gap in agreement with experiment |^ but it may not have the same amplitude. A possible 
phenomenological form of spectral function A(k, u) is A+(k, co')+A„(k, u), where A±(k, u) = 



F (1 ± ^k/^k)/ [{uJ T E^f + F^J . Here = V^k + ^k and the factors (1 ± ^k/^k) do not 
imply the phase coherence of BCS theory. 

If the frequency cut-off ujc is much larger than any energy scale in our consideration, 
then under the assumption of a cylindrical Fermi surface, it can be shown that as T ^ 0, 
-T E.,k ^(k, uf = 1/2 + K (^A(0)/F,) /vr and T E.,k i^(k, ujf = 1/2 - K (zA(0)/F,) /vr 
for superconducting state with an in-plane scattering rate F^, which is assumed, for simplic- 
ity, to be frequency independent, and — T X]t<;,k G'o(k, u;)^ = 1/2 + K {iK/V^ jn for normal 
pseudogap state. Here K is the complete elliptic integral of the first kind and we can set 
the density of states A^(0) at the Fermi level equal to be unity because it does not appear in 
NMSW. The in-plane angle dependence of t^. such as cos'^(20) does not change the results 



as long as ujc ^ ^(0) and A. Therefore, NMSW becomes 

(iV. - N.) 1 , K(#)-K(i^) 

P. 2+ ^_2K(iMa) ■ 

may be different from F and much smaller as ARPES data imply 0. Such data are 
more consistent with the existence of quasiparticles in the superconducting state than in the 
normal state. 

It is worthwhile illustrating the implications of Eq. @) in some detail, i) Suppose the 
in-plane scattering rate Tg is negligible compared with A(0), then Eq. (^) reduces to {Nn — 
Ns)/Ps = 1/2 + K (iA/T^ /tt, which depends only on the normal state and NMSW > 1/2. In 
this case, the conventional sum rule is recovered when A/F ^ 1 because the spectral function 
A(k,uj) becomes insensitive to the pseudogap. ii) If A(0)/Fs ^ 1 and A/F ^ 1, then 
NMSW ~ 1/2. To get 1/2, it is not necessary that A(0) ~ A and F^ ~ F simultaneously; in 
other words, a mismatch between A and A(0) does not matter as long as F and F^ are both 
negligible. But a small value of F gives a frequency dependence of a'^^{uj), which disagrees 
with the observation that it is flat for a; up to a few 100 meV. Thus, the large F and F^ 
limit must be applied, and consequently any pseudogap region at low uj ^ A will be filled in 
and the model cannot describe the observation made in underdoped YBCO. A cancelation 
between Z]<^,k G'(k, c<j)^ and Z]<^,k G'o(k, u;)^ in Eq. (H) can still arise for large F if there is 
a match between A(0) and A and between F and F^ as can be seen from Fig.l, where we 
plot — T X]t<;,k G'o(k, cu)^ (or —TJ2uj,]iG(k,uy) for three values of the pseudogap (or gap) as 
a function of F (or F^). We also show T X^Lj.k -^(k, cu)^. Clearly the large F region can also 
give a value of the sum rule bigger or smaller than 1/2. Should the pseudogap be larger 
than the superconducting gap (point b in Fig. 1) and F^ not too much smaller than F (point 
a), the sum rule will be less than 1/2 while if F^ is much less than F such as for point d, 
it will be larger than 1/2 but less than 1. Other combinations could also be used. To be 
more explicit, suppose that A/F = A(0)/Fs + 7 with A(0)/Fs ^ 7, then K(zA/F)/7r — 
K(^A(0)/F,)/7r ^ -(7/4)A(0)/F,. Since 1/2 - K(«A(0)/F,)/7r ^ (1/8)(A(0)/F,)2 for 
A(0)/F, < 1, NMSW^ 1/2 - 7F,/A(0) for 7 < A(0)/F, < 1. Based on ARPES data, we 
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may assume that A =i A(0) and T > ^ 0. In this case NMSW~ 1/2 + (a - d)/(2c), 
which is obviously greater than 1/2, but small Tg will not give the observed flat response 
crlci^) in the superconducting state. 

So far we have assumed that both superconducting and normal state are at T ^ 0. 
Because it is not possible to access the normal state at low T, Basov et al. have used T = Tc 
instead. In the inset of Fig. 1, we plot — T'X^cjk^o with A = 4Tc and F = Tc as a function 
of T. (We emphasize that -T E<^,k Gl is nothing but NMSW if F, < A(0).) It is clear that 
including the T dependence in the NMWS does not qualitatively change the physics of the 
conductivity sum rule. 

Incoherent coupling. For incoherent c-axis coupling, we need a specific model for 
the impurity scattering potential |Vk-pP and need to average over impurity configura- 
tion. Part of the disorder scattering can be due to a mismatch of overlap matrix ele- 
ments t± between planes. We use a simple model for the scattering potential |Vk_pp = 
+ l^iP cos(20fc) cos(2(/)p) and assume = |V^ip for simphcity. The NMSW is now 

Ps 2 2E^[(«:'V«:)K(K)-E(K)/«:f' 

where k = A/^A^ + (cu + Fsgncj)2, k' = VT^^, k = A{T)/^A{Tf + [uo + F,sgna;)2, 
k' = \J\ — K^^i and E is the complete elliptic integral of the second kind. As one can easily 
see, 'k!'^Y^{K) becomes [t^ jl^' if A = regardless of F. It can also be inferred from Eq. @ 
that i) as A increases or F decreases, NMSW becomes smaller because k is closer to 1 for 
a given tu, and ii) as F^ increases NMSW becomes smaller because n is closer to 0. In the 
incoherent coupling case, the normal state response 0""^ ^ I^P is independent of uj and there 
is no restriction on F @J^. In the inset of Fig. 2, we plot NMSW as a function of F^ with 
fixed values of a = and 1, where a = A/A(0), and F = A(0). As expected, NMSW with 
a = 1 decreases with increasing F^. Note that NMSW with a = is almost unchanged with 
increasing F^ because the numerator and the denominator of the second term in Eq. are 
both decreasing in similar manner. In the case A = A(0)(1 + 5) \b <^ 1] with F = F^ = 0, 
Eq. (D reduces to NMSW^ 1/2 - 1.083 5, which is less than 1/2. For F = F, ^ 0, the 
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coefficient of S will be changed. In the general case we need to evaluate the sums in Eq. |^ 
for given A, A(0), F and Tg. In the main frame of Fig. 2, we plot NMSW as a function of 
F with a few values of a assuming F<j = since we have seen that F^ 7^ decreases NMSW. 
As one sees, NMSW is greater than 1 and independent of F if a = 0. Note that NMSW 
with a = is an asymptote of NMSW with nonzero a because the effect of A on NMSW 
vanishes when F ^ A. As we expect, NMSW decreases with increasing a and decreasing F. 
If a > 1.46 and F ^ A(0), then NMSW< 0. Generally speaking, in the incoherent coupling 
model, we can explain NMSW^ 1/2 with A ^ A(0), F > A(0) and F, < A(0). 

Additional insight into the sum rule can be obtained from consideration of aij:'^\u!). In 
the top panel of Fig. 3, we compare (solid curve), with pseudogap A = A(0), with 

cT^^(co') (dashed curve). In both cases, F is taken to be zero and therefore both curves go 
to zero at uj = 0. At small u ^ A(0) the solid and dashed curves are almost the same. 
However, they start to differ at a higher u, with crlci^) falling below ai^^u). We already 
know that for these parameter the sum rule is exactly 1/2. But in this example, the effective 
NMSW defined as Af{uj) = [Nn{uj) — Ns{u)] / ps is nealy zero up to > A(0) as shown in 
the dashed curve in the bottom panel of Fig. 3. This does not agree with the findings of 
Basov et ai, where //{u) increases almost linearly out of zero and rapidly move toward its 
saturated value. It is clear that, even for the incoherent c-axis coupling case, the argument 



presented in Ref. |jT3[ that the preformed pair model agrees with a sum rule of 1/2 due to 
a simple cancelation between the normal and superconducting state is deficient when its 
detailed approach toward its saturated value is considered, i.e. when we compare the low lj 
behavior of ^/{uj) with experiment. 

A more reasonable model is obtained for finite normal state F. Returning to the solid 
curve of the top panel and including a finite F in the pseudogap state fills in the region 
below 2A by transfering spectral weight out of the region above this. In the middle panel 
we show results (solid curve) for a pseudogap A = 1.4A(0) and F = 0.5A(0) for illustration 
only. Other choices could also have been made showing that our explanations are robust. 
Comparison of the sohd curve (normal state) with the superconducting (dashed) curve (F^ = 



0) shows a very different low uj behavior. Now M{uj) will grow rapidly as uo increases from 
uj = 0. This behavior is shown as the solid curve of the bottom panel. M{uj) grows linearly 
out of cij = and nearly reaches its saturated value(^ 0.52) by ^ 6 ~ 7A(0) in much 
better agreement with experiment. One can also conceive that A and F are larger than 
considered above and that the normal state may not exhibit much of a pseudogap because 
of a large F, but that the sum rule is still nevertheless less than 1. This could explain 
optimally doped T12201 with NMSW~ 0.6 where, however, no in-plane pseudogap has been 
reported although in the data of Basov et al, cri^{uj) is droping with decreasing u at small 
uj consistent with a smeared pseudogap. 

Conclusion. We conclude that incoherent c-axis coupling can more easily explain the 
observed violations of the optical sum rule than can coherent coupling, but that non-Fermi 
liquid in-plane behavior is required. It is possible that the normalized missing spectral 
weight is less than 1/2. A sharp increase in accumulated normalized spectral weight as a 
function of cut-off frequency rising to a value close to its saturated value within a few times 
the gap as is observed can also be understood within our phenomenological model but not 
within a pure preformed pair model 
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FIGURES 

FIG. 1. The quantity -T Y.^ -^^ G§ (or -T Y.^^-^^ G^) as a function of T/Tc (or r^/Tc) for different 

values of A and A(0): A = Tc (dashed curve), A = A(0) = 2.5Tc (solid curve) and A = 4Tc 

(dotted curve). The solid curve is TYiui,)i^^- The inset shows the T dependence of —TYiu>,kGo 

with A = 4Tc and T = T^. 

FIG. 2. NMSW as a function of r/A(0) for different values of a = A/A(0) with = and 

T ^ A(0). The inset gives NMSW as a function of F5/A(0) for two specific cases: a = with any 

r (dashed curve) and a = 1 with F = A(0) (solid curve). 

FIG. 3. Top panel shows ai^{uj)/acn with A = A(0) and F = (solid curve), where 

(Ten = 47mj (e(iiV(O)Vb)^ with an impurity density n^. The dashed curve {ai^{Lo)/acn is for the super- 
conducting state with A(0) and F^ = 0. Middle panel gives a'^^{ijj)/acn but now with A = 1.4A(0) 
and F = 0.5A(0). Note the transfer of spectral weight from higher to lower frequency as compared 
to case when F = (solid curve in the top panel). Bottom panel shows the effective NMSW N{u). 
The dashed curve is for the preform pair model with F = F^ = and A = A(0). Its saturated 
value of 1/2 is not reached at a; = 6A(0) where it is still only about 0.3. The solid curve includes 
smearing into the normal state with A = 1.4A(0) and F = 0.5A(0). Its saturated value is about 
0.52 which is nearly reached by a; = 6A(0). 



10 



Fig. 1 (Kim & Carbotte) 
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Fig. 2 (Kim & Carbotte) 




Fig. 3 (Kim & Carbotte) 
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